% AssemblageRigidite retourne la matrice de rigidité globale


function K = AssemblageRigidite(maillage,poutre)

K = zeros(2*length(maillage));

for i=1:(length(maillage)-1)
    
iCourant = (i-1)*2+1;
elementCourant = maillage(i:i+1);
    
K( iCourant:iCourant+3 ,iCourant:iCourant+3 ) = K( iCourant:iCourant+3,iCourant:iCourant+3 ) + RigiditeElementaire(elementCourant,poutre);

end
 
end
